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ABSTRACT 

We calculate non-axisymmetric oscillations of uniformly rotating polytropes magne¬ 
tized with a purely toroidal magnetic field, taking account of the effects of the de¬ 
formation due to the magnetic field. As for rotation, we consider only the effects of 
Coriolis force on the oscillation modes, ignoring those of the centrifugal force, that is, 
of the rotational deformation of the star. Since separation of variables is not possi¬ 
ble for the oscillation of rotating magnetized stars, we employ finite series expansions 
for the perturbations using spherical harmonic functions. We calculate magnetically 
modified normal modes such as g-, /-, p-, r-, and inertial modes. In the lowest or¬ 
der, the frequency shifts produced by the magnetic field scale with the square of the 
characteristic Alfven frequency. As a measure of the effects of the magnetic field, we 
calculate the proportionality constant for the frequency shifts for various oscillation 
modes. We find that the effects of the deformation are significant for high frequency 
modes such as /- and p-modes but unimportant for low frequency modes such as g-, 
r-, and inertial modes. 
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1 INTRODUCTION 

Since the first discovery report of quasi-periodic oscillations (QPOs) in the tail of the giant X/ 7 -ray flare from the soft 7 -ray 
repeater (SGR) 1806-20 (Israel et al. 2005), extensive theoretical studies have been carried out to identify physical mechanisms 
responsible for the QPOs. SGRs belong to what we call magnetars, neutron stars possessing an extremely strong magnetic 
field as strong as 10^®G at the surface. Giant flares have so far been observed only from three SGRs, that is, SGR 0526-66, 
1900-1-44, and 1806-20, and just once for each of the SGRs, indicating that giant flares from magnetars are quite rare events. 
For magnetars, see reviews, for example, by Woods & Thompson (2006) and Mereghetti (2008). Analyzing archival data of 
another magnetar candidate SGR 1900-1-14, Strohmayer & Watts (2005) have succeeded in identifying QPOs in the X-ray 
giant flare that was observed in 1998. QPOs frequencies now identified in the giant flares from the two magnetar candidates 
are 18, 30, 92.5, 150, 626 Hz for SGR 1806-20 (Israel et al. 2005; Strohmayer & Watts 2006; Watts & Strohmayer 2006), and 
28, 53.5, 84, 155 Hz for SGR 1900-1-14 (Strohmayer & Watts 2005). Employing Bayesian statistics, Hambaryan et al. (2011) 
have reanalyzed the data for SGR 1806-20 to identify QPO frequencies at 16.9, 21.4, 36.8, 59.0, 61.3, and 116.3 Hz. For the 
giant flare in 1979 from SGR 0526-66, Watts (2011) mentioned in her review paper a report of a QPO at ~ 43Hz, but she also 
suggested the difficulty in the frequency analysis in the impulsive phase of the burst. Here, it is worth mentioning a promising 
recent attempt to find QPOs in short recurrent bursts in SGRs. Huppenkothen et al. (2014) have succeeded in identifying 
candidate signals at 260, 93, and 127 Hz from J1550-5418, where they used Bayesian statistics for the analysis. 

The QPOs are now regarded as a manifestation of global oscillations of the underlying neutron stars, and it is expected 
that they can be used for seismological studies of the magnetars. Seismological studies of magnetars may have started with a 
paper by Duncan (1998), who suggested that frequent starquakes in SGRs could excite crustal toroidal modes and the burst 
emission modulated at the mode frequencies would be detected. The detection of QPOs in the giant flare from SGR 1806-20 in 
2004 (Israel et al. 2005) has given a huge trigger leading to subsequent intensive theoretical studies of QPOs in magnetars. In 
the early studies of QPOs in magnetars (e.g., Piro 2005; Lee 2007), the oscillations were assumed to be practically confined in 
the solid crust, as Duncan (1998) anticipated, and the effects of a magnetic field in the fluid core were ignored. Since magnetars 
are believed to possess an extremely strong magnetic field threading both the solid crust and the fluid core, to determine 
the oscillation frequency spectra of the stars, we need to correctly take account of the effects of the strong magnetic field in 
both regions on the oscillations. Applying a toy model, Glampedakis, Samuelsson, and Andersson (2006) discussed toroidal 
oscillations as global discrete modes residing in the fluid core and in the solid crust both threaded by a strong magnetic field, 
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and showed that the modes that are most likely to be excited in magnetars are such that the crust and the core oscillate 
in concert. Also using a toy model, Levin (2006, 2007) put forward a different idea that Alfven modes in the fluid core may 
lead to the formation of continuum frequency spectra and that toroidal crust modes will be rapidly damped as a result of 
resonant absorption in the core, and even suggested that there exist no discrete normal modes for the strongly magnetized 
neutron stars. However, assuming a pure poloidal field threading both the solid crust and the fluid core, Lee (2008) and Asai 
& Lee (2014) carried out normal mode calculations of axisymmetric toroidal modes and found discrete toroidal modes. Later 
on, van Hoven & Levin (2011, 2012), employing spectral method of calculation, have succeeded in suggesting the existence 
of discrete modes in the gaps between continuum frequency bands. In general relativistic frame work, Sotani et al. (2007) 
computed normal modes of magnetized neutron stars in the weak magnetic field limit. 

Oscillations in magnetized stars are governed by a set of linearized partial differential equations. Normal mode analysis 
of the oscillations of magnetized stars are not necessarily easy to conduct, partly because separation of variables between the 
radial and the angular coordinates is in general impossible to represent the perturbations, and partly because the possible 
existence of continuum bands in the frequency spectra make it difficult to properly calculate oscillation modes, particularly 
belonging to the continua (see, e.g., Goedbloed & Poedts 2004). In normal mode calculations, we usually employ series 
expansions of a finite length to represent the perturbations so that the set of linearized partial differential equations reduces 
to a set of linear ordinary differential equations for the expansion coefficients. This process can be very cumbersome when 
we try to carry out modal analyses for various configurations of magnetic fields. This may be one of the reasons why MHD 
simulations have been used to investigate the small amplitude oscillations of magnetized neutron stars by many authors 
including, e.g., Sotani, Kokkotas & Stergioulas (2008), Cerda-Duran et al. (2009), Colaiuda & Kokkotas (2011, 2012), Gabler 
et al. (2011, 2012, 2013a,b), Lander et al. (2010), Passamonti & Lander (2013, 2014). In the analyses with MHD simulations, 
QPOs are believed to be associated with continuum spectra and should be properly distinguished from discrete normal modes 
by closely watching motions and phases of various points in the interior. 

Gonfigurations of magnetic fields in magnetars are highly uncertain (e.g., Thompson & Duncan 1993, 1996; Thompson, 
Lyuitikov & Kulkarni 2002). As shown by the core-collapse supernova MHD simulations (e.g., Kotake, Sato & Takahashi 
2006), toroidal fields can be easily produced and amplified when winding of the initial seed poloidal fields is effective in the 
differentially rotating core even if there is no initial toroidal fields. For modal analyses, it is desirable to examine various 
magnetic field configurations. Most of the modal analyses so far carried out have been for a purely poloidal magnetic field 
(Lee 2007, 2008; Sotani et al 2007, 2008; Gerda-Duran et al. 2009; Golaiuda & Kokkotas 2011, 2012; Gabler et al. 2011, 
2012; Passamonti & Lander 2013, 2014; Asai & Lee 2014). Recently, however, some authors, using MHD simulations, started 
investigating small amplitude oscillations for a purely toroidal magnetic field (Lander et al. 2010; Passamonti & Lander 2013), 
and even for mixed poloidal and toroidal field configurations (Gabler et al. 2013). Using MHD simulations, for example. 
Lander et al (2010), calculated rotational modes (r-modes and inertial modes) of magnetized stars and showed that r-modes 
at rapid rotation tend to magnetic modes in the slow rotation limit. 

In this paper we carry out normal mode analysis of polytropic models with a purely toroidal magnetic field for various 
non-axisymmetri oscillation modes, including p-, /-, p-modes, and rotational modes such as r-modes and inertial modes, where 
we include the effects of equilibrium deformation due to the magnetic field on the oscillations. To calculate rotational modes, 
we only take account of Coriolis force and include no effects of the centrifugal force. The oscillation equations for magnetically 
deformed rotating stars are derived by following the formulation similar to that by Saio (1981) (see also Lee 1993; Yoshida 
& Lee 2000a). The numerical method to compute normal modes of magnetized rotating stars is the same as that in Lee 
(2005) (see also Lee 2007), who employed series expansions of a finite length in terms of spherical harmonic functions for 
the perturbations. This paper is organized as follows. §2 describes the method used to construct a magnetically deformed 
equilibrium stellar model, and perturbation equations for non-axisymmetric oscillation modes in magnetized rotating stars 
are derived in §3. Numeical results are summarized in §4 and we conclude in §5. The details of the oscillation equations solved 
in this paper and suitable boundary conditions imposed at the stellar center and surface are given in Appendix A. 


2 EQUILIBRIUM MODEL 

We consider the oscillations of uniformly rotating polytropes with purely toroidal magnetic fields. Equilibrium structures of 
stars having purely toroidal magnetic fields have so far been studied with non-perturbative approaches within the framework 
of Newtonian mechanics (Miketinac 1973) and of general relativity (Kiuchi & Yoshida 2008; Frieben & Rezzolla 2012). In 
this study, we employ a perturbative approach to construct stars deformed by a purely toroidal magnetic field. Following 
Miketinac (1973), a purely toroidal magnetic field imposed on the stars in equilibrium is assumed to be given by 

Br — Q, Be=Q, B,/, = kprsinO, (1) 

where k = Bo/{y/ 2 pcR) is a constant. Bo is the parameter used for the strength of the magnetic field in the interior, pc 
is the density at the stellar center, and R is the radius of the star. Here, we use spherical polar coordinates {r, 6 ,(j>). The 
magnitude of the magnetic field is given by |B| = \BA = {Bo/V2)pxsm6, where x = r/R and p = pjpc- The fluid velocity v 
in equilibrium is assumed to be given by 


V-. 


0 , ve = 0 , V 4 , = rsinfffl, 


( 2 ) 
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where Q. denotes the angular velocity of the uniformly rotating star. In this study, the deformation of the star is assumed to 
be solely caused by the magnetic fields and the effects of the centrifugal force are ignored. By the assumptions (1) and (2), 
the induction and continuity equations are automatically satisfied and need not be considered any further. For the toroidal 
field (1), we can write the Lorentz force per unit mass as a potential force, that is, 

X B) X B = -V (^-px^ sm . (3) 

47rp \ Svrpc J 

The structure of a star in equilibrium is then determined by the hydrostatic equation, the Poisson equation, and the equation 
of state: 


Vp = —(4) 
= drrGp, (5) 


P = 


it'cP 


l + l/n 


( 6 ) 


where n and Kc are the polytropic index and the structure constant given by the mass and the radius of the star, G is the 
gravitational constant, 'F is the gravitational potential, and is the effective potential defined by 

4'= $ + icuir^p [1 - P 2 (cos 6 »)] - G, (7) 

where wa = \//[AirpcR^) is the characteristic Alfven frequency of the star and G is a constant. Here, P 2 (cos 6 ) = (3 cos^ 6 — 
l)/2 denotes the Legendre polynominal of order 2. 

Since the potential 4> is the quantity of order of GM/R with M being the mass of the star, the ratio of the second term 
on the right hand side of equation m to $ may be given by oj\ = where Hx = •JGMjR?. For a neutron star model 

of the mass M — IAMq and radius R — 10®cm, for example, we have uj\ ~ 2 x 10“® for the field strength Bo = 10^®G, 
suggesting that the effects of the magnetic held on the equilibrium structure is not signihcant so long as Bo < lO^^G. In this 
paper, as mentioned before, we assume that the magnetic held is sufhciently weak so that the deformation of the equilibrium 
structure due to the magnetic held can be treated as a small perturbation to the non-magnetic stars when ijj\ 1. Under this 
assumption, we can regard p appearing in the terms proportional to in equation 0 as the density po in the non-magnetic 
star. Thus, T satishes 


= 47rGp-f itui 


2 f>o 


+6r^+6po- 




Since p can be regarded as a function of from equations 0 and ®, if we expand T (r, 6 ) as 
\['(r, 6 ) = To(r) - 2 R^uj\ [^o{x) + ■i/' 2 (x)P 2 (cos 0 )] , 


( 8 ) 

( 9 ) 


we may expand p(r, 6 ) as 

p{r,e) = po(r) - [i>o{x) -|- V’ 2 (a:)P 2 (cos 6 »)]. ( 10 ) 

QWo 

Here, 4'o(r) = 4?o(r) and po(r) are the gravitational potential and the density of the non-magnetized star, and they satisfy 
dpo/dr = —pod^o/dr = —poGMr/r'^, Mr = dnr^podr, and po = 

Substituting equations and GQl) into 0, we hnd 

R^V^ [V’o(*) + ip2(x)P2(cos0)] = 47 rGP^^^ bPo{x) -|- 'ip 2 {x)P 2 {cosd)] + fo{x) + f 2 {x)P 2 {cosd), ( 11 ) 

from which we obtain the following linear ordinary differential equations for ipo{x) and 'tp 2 {x): 


_L A 

x^ da: \ da; / 


k{x) + 


i>2 + f2ix), 


(13) 


where 


foix) 


1 

6 



-I- -I- 6po 

dr 


h{x) 


1 

6 



-I- 6r 


dpo 

dr 


(14) 


and 


k{x) = inGR^^. 


(15) 
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In order to numerically integrate the differential equations m and m from the stellar center, we need to impose the 
regularity condition at the center, which may be obtained by substituting the expansion around the center x = 0 


ipi = *" 


(16) 


into m and m for j — 0 and j = 2, respectively. Since k{x) —>■ fc(0), foi^:) /o(0), and f 2 {x) f 2 oY as ® 0, where 

A:(0), /o(0), and /20 are constants, values of the exponent s are given by s = j for the regular solution at the center and the 
expansion coefficients and remain undetermined. Assuming that the density at the center x = 0 is independent of u)a, 
we have ag°^ = 0 for 'ipQ. The expansion coefficient for 'ip 2 must be specified by applying the surface boundary condition: 


See Appendix B for the derivation of the boundary condition. 


(17) 


3 PERTURBATION EQUATIONS 

For the modal analysis of magnetically deformed stars, we introduce the parameter a that labels equi-potential surfaces of 
^(r, 0). The parameter a is defined such that ’^(r,9) = ’Fo(n), that is, 


To(a) = 'l'o(r) - 2 R^ [V’o(*) + V’2(x)P2(cos6»)], 


(18) 


which may define the equipotential surface as given by a function r(a, 9). Assuming the deviation of the equipotential surface 
r = r(a, 9) from the spherical surface r = a is small, we dehne the function r(a, 9) as 


r = a [1 + e(a, 6>)] , 


(19) 


and we assume that e is the quantity of order of By substituting equation (I19II into (II8II . we obtain the explicit 

expression for the function e{a,9) up to the order of 

e{a,9) = a{a) + /3{a)P2{cos9), (20) 

where 

a{a) = /3(a) = V’2(a:), (21) 

where ci = {a/RA/[M{a)/M], and M{a) denotes the mass inside the a-constant surface and M = M{R). 

Hereafter, we employ the parameter a instead of the polar radial coordinate r as the radial coordinate. In this coordinate 
system (a, 9, <j)), the line element is given by 

ds^ = (1 + 2e)(da^ + d9^ + sin^ 9d(jL) + 2a'^da^ + 2a'^dadS. (22) 

Od Ou 

Note that in this coordinate system, the pressure, the density and the effective potential of a magnetized star depend only on 
the radial coordinate a, although the orthogonality of the basis vectors is lost. 

The governing equations of non-radial oscillations of a magnetized and uniformly rotating star are obtained by linearizing 
the basic equations. As for rotation effects on the oscillations, as mentioned in the previous section, we consider only the effects 
of the Coriolis force and ignore those of the centrifugal force, where we assume the rotation axis is parallel to the magnetic 
axis. Since the equilibrium state is assumed to be stationary and axisymmetric, the perturbation quantities are proportional 
to exp(ia;t + im(f>), where w is the frequency observed in an inertial frame and m is the azimuthal wave number. Then, the 
linearized basic equations which govern the adiabatic, non-radial oscillations of a magnetized and uniformly rotating star are 
written in the coordinate system {a,9,(j)), to second order in loa, as (Saio 1981; Lee 1993; Yoshida & Lee 2000a) 


- (7^ [(1 -b 2£)^ -b ae“Voe -b a(^ • Voe)e,] = -Vo$' --Vop + 


^{Vox B) X B 
47r 


da 


-b iaD 


+ — [(Vo X S') X -B -b (Vo xB)x B'] , 

47rp '■ 


p' + Vo • (p^) + p^ ■ Vo (^3e + a|^^ = 0, 


(23) 


(24) 


P_ 

P 


f - ^-.A. 

1 IP a 


(BT = 


ijk 


d 


y/g dx^ 




(25) 

(26) 
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where a = tj + mO. denotes the oscillation frequency observed in the corotating frame of the star, ^(a, 6 , (j>) is the displacement 
vector, the prime {') indicates the Eulerian perturbation, eijk and are the Levi-Civita permutation symbols, g is the 
determinant of the metric gij, and 


Vo = lim 

e-fO 


d Id Id 

e-a -1" ^—- — 

da a da a sin o dtp 


(27) 


and Ga, 60 , and are the basis vectors in the a, 0, and (j) directions, respectively. Here, the vector D in equation (|23p conies 
from the Coriolis force and is given by (see, e.g., Lee 1993; Yoshida Sz Lee 2000a) 

Da = 2Q(l + 2e + a^) sin , Dg = 2n(l + 2e+ ^ cos , 

\ da) \ cos 6 dO 


Ds = — 2 fl 


1 r. de\ , a / sin0 de 

1 + 2e + a^— ) sm 9^ + I 1 + 2e H-- — 

da J V cos 0 dO 


cosOC 


and aA in equation (I25II denotes the Schwarzschild discriminant defined as 
dlnp 1 dlnp 


aA = 


d In a r 1 d In a ’ 


(28) 


(29) 


where Fi = (91np/91np)ad- In this paper, for simplicity, we employ the Cowling approximation, neglecting il?'. 

Because of the Lorentz and Coriolis terms in the equation of motion (1231) . separation of variables for the perturbations 
is impossible between the radial coordinate (a) and the angular coordinates {9, <j>). We therefore expand the perturbations in 
terms of the spherical harmonic functions Y{"{9, (j>) with different Vs for a given azimuthal index m. The displacement vector 
^ is then given by (see e.g., Lee 2005, 2007) 




(30) 


j = l L 


-a, (“) ® V)".« - “r,; (9. ♦) 


(31) 


^^ = E 

J =1 




and the vector B' is given by 

/Day 

= E iahffia)Yi^{9,(l)), 


(32) 


(33) 


(B 


kp 


0\/ Jmax r 

= E 

j = l L 


(a) Ay-(e, ^) _ (a)^ Ay-(g, 


(34) 


/ 1D'P\^ Jmax 

kp ^ 


^A ahj, {a) ^Y,T{9, cj.) 


(35) 


where Ij = \m\ + 2 {j — 1 ) and !'■ = Zj + 1 for even modes, and Ij = \m\ + 2ji — 1 and Z) = Zj — 1 for odd modes, respectively, 
and j = 1, 2, 3,..., jmax- 

The Euler perturbations of the pressure and density are given by 


P' P'h = E P'h (“)^i 7 (^’ '(')• 

3=1 


(36) 


In this paper, we usually use jmax = 12 to obtain solutions with sufficiently high-angular resolution. Substituting the expansions 
(O- dMI into the linearized basic equations (1231) - (1261) . we obtain a finite set of coupled linear ordinary differential equations 
for the expansion coefficients Si^ (a) and p[. (a), which we call the oscillation equations to be solved in the interior of magnetized 
and uniformly rotating stars. The set of oscillation equations obtained for the magnetized rotating star is given in Appendix A. 
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Figure 1. Equi-magnetic field strength contours (solid lines) and equi-density contours (long-dashed lines) on the meridional cross 
sections are plotted, from left to right panels, for the polytropes of the index n =1, 1.5, and 3, respectively. The outer-most solid circles 
show stellar surfaces. The solid contours correspond to fJ/Bmax = 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, and 0.9, and the long-dashed 
contours to p/pc = 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, and 0.9. 



Figure 2. Functions a' = and /?' = 


versus the fractional radius a/R for the n = 1 polytrope. 


4 NUMERICAL RESULTS 

In this study, the three polytropes with indices n = 1, 1.5, and 3 are used for modal analyses of the magnetized stars. The 
polytropes with n = 1 and 1.5 and with n = 3 are regarded as simplified models of the neutron star and the normal star, 
respectively. For these polytopes, distributions of the density and the imposed magnetic fields are shown in Fig [T] In this 
figure, the equi-density and equi-magnetic field strength contours on the meridional cross sections are given. We see that the 
density and magnetic field distribution of the models with a larger polytropic index tend to be concentrated in the central 
region of the star. Typical values of the mass M and radius R for the neutron star and the normal star are assumed to be 
{M,R) — (1.4M©, 10®cm) and {M,R) = (M©,!?©), respectively. Thus, we have uia = 4.42 x 10“^ for the neutron star with 
the field strength Bo = 10^®G, and oja ~ 7.39 x 10”“^ for the normal star with the field strength Bq = 10®G. For the polytrope 
with n = 1, the functions a.(a)lY\ and P{a)lijjA for the magnetic deformation are plotted as functions of a/R in Fig. [2] 


4.1 g-, /-, and p-modes 


We first calculate /-, low radial order g-, and p-modes of the polytropes taking account of the effects of the toroidal magnetic 
field. In these calculations, no effects of rotation are considered. To study oscillation modes for the neutron star and normal 
star models, the adiabatic indices for perturbations are assumed to be 


1 

IT 


n 

n -b 1 


+ 7 


(37) 


with 7 being a constant, for which aA = 7 (cilnp/cilna). In this subsection, we use 7 = —10“"* for the polytropes with the 
indices n = 1 and 1.5, but for the polytrope with n = 3, we assume Fi = 5/3, and hence 7 = —3/20. Since all the magnetic 
terms in the oscillation equations are proportional to cLa, an oscillation frequency of the modes may be written by (see 
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Table 1. CoefBcients E 2 , E^, and E^ for g-, /-, and p-modes ot I = m for the polytropic model with n = 1 and 7 = —10 ^ * 


mode UQ E 2 E 2 i ?2 


m = 1 


93 . 0.00570 -1.578(+0) -1.490(+0) -1.491(+0) 

92 . 0.00770 -2.051(+0) -2.095(+0) -2.095(+0) 

91 . 0.01203 -3.003(+0) -3.063(+0) -3.064(+0) 

pi . 3.26931 6.162(-1) 6.164(-1) 3.553(-l) 

P2 . 5.09325 1.113(+0) 1.114(+0) 5.650(-l) 

P3 . 6.85013 1.550(+0) 1.554(+0) 7.586(-l) 


m = 2 


93 . 0.00884 3.518(+1) 3.488(+l) 3.487(+l) 

92 . 0.01152 2.427(+l) 2.402(+l) 2.402(+l) 

91 . 0.01678 1.133(+1) 1.130(+1) 1.130(+1) 

/ . 1.65562 4.107(-1) 4.109(-1) 3.567(-l) 

pi. 3.79225 7.099(-l) 7.103(-1) 4.387(-l) 

P2 . 5.67886 1.173(+0) 1.175(+0) 6.713(-1) 

P3 . 7.48089 1.605(+0) 1.608(+0) 8.874(-l) 


m = 3 


93 . 0.01136 7.583(+l) 7.498(+l) 7.498(+l) 

92 . 0.01443 5.352(+l) 5.319(+1) 5.319(+1) 

91 . 0.02008 2.732(+l) 2.732(+l) 2.731(+1) 

/ . 1.97094 6.696(-l) 6.696(-l) 5.545(-l) 

pi. 4.22956 8.280(-l) 8.282(-l) 5.196(-1) 

P2 . 6.18783 1.272(+0) 1.272(+0) 7.593(-l) 

P3 . 8.04091 1.701(+0) 1.703(+0) 9.900(-l) 


* We use the notation of 1.000 X 10^ = 1.000(77). 


Appendices A & C and Unno et al. 1989) 

d = d-Q + E 20 JA H-> (38) 

where do is the oscillation frequency of the non-magnetized star, and E 2 is a proportionality coefficient and can be obtained 
by calculating the oscillation frequency of the mode for two different values of d)^, that is, = 0 and ~ 10 ~®, for example. 
Here, do and E 2 are the quantities normalized by the Kepler frequency of the star flfc. This coefficient E 2 for a mode may 
also be calculated by using the eigenfunctions of the non-magnetized star by treating as a small parameter. We have used 
the symbol E 2 to denote the coefficient computed by using the eigenfunctions for the non-magnetized star. The derivation 
and the explicit expression for E 2 are given in Appendix C. 

In Tables 1-3, we tabulate the coefficients E2 and E2 as well as the frequency do for the /-modes and low radial order p- 
and 9 -modes of / = m = 1, 2, and 3 for the polytropes with the indices n = 1, 1.5, and 3. We observe that the two coefficients 
E 2 and E 2 are in good agreement with each other, except for a few very low frequency 9 -modes. In these tables, we also 
tabulate E 2 , which is the same as the coefficient E 2 but calculated ignoring all the equilibrium deformation effects due to 
magnetic stress. We see that the frequencies of the /- and p-modes are strongly affected by the equilibrium deformation, but 
the deformation effects are not very important for the 9 -modes. This property of the frequency responses to the deformation 
due to magnetic field is quite similar to that found for the rotational deformation (Saio 1981). We note that the frequency we 
obtain for / mode is consistent with that by Lander et al. (2010) for < 0.1 (because they consider the effects of the 

second order of Q). 

From Tables 1 and 2, we see modal properties of /-, low radial order 9 -, and p-modes for the neutron star models. Because 
of the small value of I7I, the frequencies do of the 9 -modes are quite low, which may be consistent with almost isentropic 
stratifications expected in the deep interior of cold neutron stars. We find the ratio E 2 l(Jo for the 9 -modes is much larger than 
the ratio for the /- and p-modes, suggesting the low frequency 9 -modes are more susceptible to the magnetic field, reflecting 
their very low frequencies of order of \/\7\- The ratio E 2 l<Jo for the 9 -modes increases with m, while the ratio for /- and 
p-modes only weakly depends on m. From Table 3, we see modal properties of /-, low radial order 9 -, and p-modes for the 
normal star model. The ratio if 2 /u'o have almost the same order of magnitudes for the 9 -, /-, and p-modes, except for the 
m = 1 9 -modes, for which the ratio is much smaller. It is also interesting to note that the magnitudes of the ratio for the 
9 -modes is of order of 0.1 (except for the m = 1 9 -modes), the value of which is much smaller than those for the 9 -modes of 
the polytropes of n = 1 and 1.5 with 7 = —10“"^. 

The coefficient E 2 in the tables 1 to 3 is computed by using two data points with different values of For example, 
E 2 computed with 4 data points for 93 modes of m = 1 and m = 2 for n = 1 are —1.577 and 35.18, and the coefficients for 
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Table 2. Coefficients E 2 , E^, and E^ for g-, /-, and p-modes oi I = m for the polytropic model with n = 1.5 and 7 = —10 ^ * 


mode CTO E 2 i ?2 


m = 1 


33 . 0.00788 -8.899(-l) -8.120(-1) -8.123(-1) 

32 . 0.01057 -1.179(+0) -1.193(+0) -1.193(+0) 

31 . 0.01626 -1.797{+0) -1.856(+0) -1.856(+0) 

pi. 3.08199 5.763(-l) 5.761(-1) 9.406(-2) 

P2 . 4.64233 9.672(-l) 9.668(-l) 1.600(-1) 

P3 . 6.15692 1.329(+0) 1.328(+0) 2.190(-1) 


m = 2 


33 . 0.01214 2.316(+1) 2.286(+l) 2.286(+l) 

32 . 0.01564 1.628(+1) 1.607(+1) 1.607(+1) 

31 . 0.02217 8.198(+0) 8.153(+0) 8.152(+0) 

/ . 1.84930 4.027(-l) 4.027(-l) 1.621(-1) 

pi. 3.55537 7.093(-l) 7.093(-l) 1.205(-1) 

P2 . 5.14850 1.073(+0) 1.073(+0) 1.808(-1) 

P3 . 6.69114 1.429(+0) 1.428(+0) 2.454(-l) 


m = 3 


33 . 0.01547 4.950(+l) 4.872(+l) 4.872(+l) 

32. 0.01935 3.546(+l) 3.511(+1) 3.511(+1) 

31 . 0.02596 1.941(+1) 1.938(+1) 1.938(+1) 

/ . 2.15084 5.640(-l) 5.640{-l) 2.301(-1) 

pi. 3.93952 8.393(-l) 8.391(-1) 1.504(-1) 

P2 . 5.58066 1.192(+0) 1.191(+0) 2.007(-l) 

P3 . 7.15896 1.544(+0) 1.542(+0) 2.653(-l) 


* We use the notation of 1.000 X 10^ = l.OOO(Af). 





Figure 3. Eigenfunctions of m = 2 even modes for the polytrope with n = 1 and 7 = —10~^ for Bq = 10^® G, where, from left to right 
panels, the eigenfunctions plotted are those of the 31 , /, and pi modes. The solid lines, the long dashed lines and the short dashed lines 
are for the functions xSi ,, xHi ,, and xTi,^i with x = ajR, and the amplitude normalization is given by Si, = 1 at the surface. 


Pi modes of m = 1 and m — 2 are 0.6164 and 0.7100. We therefore think that the coefficients E 2 in the tables have at least 
three significant digits. 

In Figure 3, we plot the expansion coefficients Si,, Hi,, and of the gi, f, and pi-modes of Z = m = 2 for the n = 1 
polytrope with Bq — 10^® G. The first expansion coefficients associated with the harmonic degree h and I'l are dominant 
over the coefficients with Ij and I'j for j > 1 , and the difference in the dominant expansion coefficients of the modes between 
the magnetized and non-magnetized models is almost indiscernible. Because of the surface boundary condition (IA14II and an 
algebraic relation (Eil) in Appendix A, when = 0 and ^ 1 we have Hi, ~ Si,/a^ at the surface and hence Hi, at the 
surface can be very large for 3 -modes having very low frequencies (f <C 1 for the normalization Si, = 1. In Figure 4, we plot 
magnetic field perturbations Sti, Hu, and Tu' of the gi, f, and pi-modes of Z = m = 2 for the n = 1 polytrope with Bq = 10 ^® 
G, where Sbi = kpahf /Bq, Hu = kpahf /Bo, and Tu' = kpahf, /Bq. 
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Table 3. CoefBcients E 2 , E^, and E^ for g-, /-, and p-modes ot I = m for the polytropic model with n = 3 and Fi = 5/3 * 


mode UQ E 2 E 2 i ?2 


m = 1 


93 . 0.88994 -1.088(-3) -7.875(-4) -7.747(-3) 

92 . 1.16154 -1.270(-3) -1.079(-3) -9.976(-3) 

91 . 1.68082 9.069(-4) 7.333(-4) -1.306(-2) 

pi . 3.81006 2.906(-l) 2.905(-l) 2.689(-3) 

P2 . 5.01208 4.335(-l) 4.331(-1) 2.803(-3) 

P3 . 6.25522 5.541(-1) 5.538(-l) 3.836(-3) 


m = 2 


93 . 1.35792 1.948(-1) 1.930(-1) 1.802(-1) 

92 . 1.70580 1.502(-1) 1.494(-1) 1.332(-1) 

91 . 2.29614 1.103(-1) 1.099(-1) 8.008(-2) 

/ . 3.06379 2.193(-1) 2.194(-1) 1.814(-2) 

pi . 4.14666 3.521(-1) 3.521(-1) 1.064(-2) 

P2 . 5.39097 4.792(-l) 4.788(-l) 7.335(-3) 

P3 . 6.65382 5.990(-l) 5.978(-l) 6.608(-3) 


m = 3 


93 . 1.70370 4.161(-1) 4.122(-1) 3.942(-l) 

92 . 2.07374 3.245(-l) 3.227(-l) 3.007(-l) 

91 . 2.64527 2.338(-l) 2.334(-l) 1.920(-1) 

/ . 3.12498 2.567(-l) 2.566(-l) 3.173(-2) 

pi. 4.37567 3.947(-l) 3.946(-l) 1.118(-2) 

P2 . 5.68481 5.193(-1) 5.189(-1) 8.675(-3) 

P3 . 6.98164 6.392(-l) 6.383(-l) 8.076(-3) 


* We use the notation of 1.000 X 10^ = 1.000(77). 



a/R 




a/R 


Figure 4. Same as Figure 3 but for the eigenfunctions = kpahf /Bq (solid lines), Hi,i^ = kpah^/Bg (long dashed lines), and 

'^ 6 * 1+1 = kpahj, /Bf) (short dashed lines). 


For slowly rotating stars, we may write the inertial frame oscillation frequency w of a mode as 

U) = CJQ “f 7Tl{Cl — l)ff + E 2 UIJI -!-•••, (39) 

where Ci represents the response of the mode frequency to the slow rotation. Since Vi\ ~ 10 “® ~ 10 “® for Bg = 10 ^® ~ lO^^G 
and E 2 ~ 10 , the rotational effects may dominate the magnetic ones for fl > 10 ”^. 


4.2 rotational modes 

We consider the effects of the magnetic field on rotational modes such as inertial modes and r-modes, for which the Coriolis 
force is the restoring force and the oscillation frequency is proportional to the rotation frequency fl. As shown by Yoshida 
& Lee (2000b), the stratification of the stellar interior strongly affects modal properties of the rotational mode. Since we are 
concerned with purely magnetic effects on the rotational mode, in this subsection, we focus on non-stratified stars or isentropic 
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Figure 5. Eigenfunctions of m = 2 rotational modes for the isentropic n = 1 polytrope for Bq = 10^® G: r mode of kq = 0.6667 (left), 
and inertial modes of kq = 1-1000 (center) and kq = 0.5173 (right). The solid lines, the long dashed lines, and the short dashed lines are 
for the functions xSi-,^, xHi,^, and with x = a/R, respectively, and the amplitude normalization is given by max(a;T(/^) = 1. 


stars, in which the adiabatic index for perturbations is given by Fi = 1 + . To represent the effects of the magnetic field on 

the rotational modes, it is convenient to use the frequency ratio k = cr/fl, where a — oj + mfl denotes the frequency observed 
in the co-rotating frame of the star, and for small values of cD^(<^ we may write 


K = Ko(fl) 


1 + T?2 


^2 


(40) 


where the coefficient kq may depend on the rotation rate and the coefficient r ^2 is a constant in the limit of —>■ 0 (see 

Appendix C). Inertial modes and r-modes of uniformly rotating isentropic polytropes were studied, for example, by Lockitch 
& Friedman (1999) and Yoshida & Lee (2000a). The r-modes of m ^ 0 and I' ^ \m\ are non-axisymmetric and retrograde 
modes and the frequency ratio kq tends to 2m/[F(Z' -f 1)] as O —>■ 0. The ratio kq for an inertial mode also tends to a definite 
value as —>■ 0, depending on m, I, and the polytropic index n (see, e.g., Yoshida & Lee 2000a), and we may use ko(0) as a 
labeling of the inertial modes for a given m. 

Since stars with strong magnetic fields are frequently very slow rotators, we consider rotational modes in slowly rotating 
stars. In Table 4, the coefficients ri 2 and 772 bs well as kq are tabulated for the I' = |m| r-modes and inertial modes for m = 2 
of isentropic (i.e., 7 = 0) polytropes with three different indices n . We use the symbol 772 to denote the coefficient computed 
by using the eigenfunctions of non-magnetized slowly rotating stars. The coefficient 772 of the fitting formula y = 772 ®, where 
X = 1 / 0 ? and y = i72/(ro, can be calculated by a least-square method. In the table, Iq — m — 1 means the r-modes, and 
lo — m ^ 2 correspond to inertial modes (see Yoshida & Lee 2000a), and the even and odd numbers of Iq — m indicate even 
and odd parity, respectively. Since using kq, the frequencies of the rotational modes can be written by k = kq -I- K 2 O?, the 
intercept kq of the fitting formula y = k,o + H 2 X can also be calculated by a least-square method. We find that the coefficients 
772 and 772 are in good agreement with each other. It is important to note that the effects of the magnetic deformation on the 
rotational modes are quite small, which is similar to the case of low frequency g-modes. We note that the frequency we obtain 
for r mode is consistent with that by Lander et al. (2010). 

In Figure 5, we show the eigenfunctions of the 7 n = 2 rotational modes of the ti = 1 polytrope for Bq = 10*^® G, where 
we assume O = 0.05 and the amplitude normalization is given by Ty^i^R) = 1. We find that the expansion coefficients for the 
m = 2 inertial mode of ko(0) = 1.1 are the same as those shown in figure 1 of Yoshida & Lee (2000a). In Figure 6, we plot 
the magnetic field perturbations Su, Hu, and Tu’ of the m = 2 rotational modes for Bq = 10^® G. In Figures 7 and 8 , we 
plot the eigenfunctions and magnetic field perturbations of the m = 2 rotational modes of the isentropic polytrope with the 
index 77 = 3 for Bq — 10® G. 


4.3 magnetic modes 

We looked for very low frequency modes to find magnetic modes having real frequencies for non-rotating stars, but we found 
none. We obtained only solutions having pure imaginary a for a given value of jmax, but we found that these solutions are 
dependent on jmax and cannot be regarded as correct solutions we look for. It is to be noted that we could not obtain very low 
frequency p-modes either in the frequency range where magnetic modes having real frequencies might coexist. The functions 
^ of the “modes” in that frequency range have discontinuities as a function of a, which suggests that the “modes” are in a 
continuum band of the frequency spectrum (see, section 7.4 of Goedbloed & Poedts 2004). This discontinuity of the functions 
which occurs in a certain low frequency range, may be caused by the relation AT' = B^ coming from equations (A4) 
and (A5) and used to eliminate the variables T” in equations (A 2 ) and (A3), where 47 = (iT,T)^ and $ = ( 2 /i, 2 / 2 )^> A 
and B are matrices, that is, there appears a point at which the determinant of the matrix A vanishes. Goedbloed & Poedts 
(2004) discussed the simplest case of a second order ordinary differential equation that possesses a frequency band in which 








magnetic perturbation 
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Table 4. Coefficient 772 of m = 2 rotational modes for isentropic polytropes with the indices n = 1, 1.5, and 3 * 


lo - I'm] 

Ko 

m 

V 2 

n = 1 

1 

0.66666 

8.324(-l) 

8.326(-l) 

2 

-0.55660 

9.105(-1) 

9.171(-1) 


1.10002 

2.398(-l) 

2.401(-1) 

3 

-1.02590 

3.304(-l) 

3.316(-1) 


0.51734 

1.784^0) 

1.785^0) 


1.35777 

1.402(-1) 

1.404(-1) 

4 

-1.27290 

2.481(-1) 

2.511(-1) 


-0.27533 

7.103^0) 

7.100^0) 


0.86296 

5.805(-l) 

5.820(-l) 


1.51956 

9.950(-2) 

9.957(-2) 

n = 1.5 

1 

0.66666 

5.247(-l) 

5.244(-l) 

2 

-0.69650 

3.688(-l) 

3.696(-l) 


1.06257 

1.752(-1) 

1.753(-1) 

3 

-1.12782 

1.439(-1) 

1.443(-1) 


0.53564 

1.071^0) 

1.069^0) 


1.31001 

1.103(-1) 

1.104(-1) 

4 

-1.34198 

8.378(-2) 

8.382(-2) 


-0.36425 

2.999^0) 

2.993^0) 


0.85864 

3.516(-1) 

3.520(-l) 


1.47217 

8.210(-2) 

8.237(-2) 

n = 3 

1 

0.66667 

1.329(-1) 

1.328(-1) 

2 

-1.07669 

2.505(-2) 

2.520(-2) 


0.99492 

3.902(-2) 

3.892(-2) 

3 

-1.37189 

1.849(-2) 

1.872(-2) 


0.57976 

2.620(-l) 

2.622(-l) 


1.20940 

2.680(-2) 

2.665(-2) 

4 

-1.51785 

6.058(-3) 

6.101(-3) 


-0.66228 

2.063(-l) 

2.064(-l) 


0.85853 

6.254(-2) 

6.253(-2) 


1.36256 

2.973(-2) 

3.023(-2) 


* We use the notation of 1.000 x 10^ = l.OOO(W). 



Figure 6. Same as Figure 5 but for the eigenfunctions Sbii = kpahf /Bq (solid lines), Hbi^ = kpahf^/B q (long dashed lines), and 
= kpahji/ Bq (short dashed lines). 
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Figure 7. Eigenfunctions of m = 2 rotational modes for the isentropic n = 3 polytrope for Bq = 10^ G: r mode of kq = 0.6667 (left), 
and inertial modes of kq = 0.9949 (center) and kq = 0.5798 (right). The solid lines, the long dashed lines, and the short dashed lines are 
for the functions xSi^, and xTi/^ with x = a/R^ respectively, and the amplitude normalization is given by max(a:T;/^) = 1. 



Figure 8. Same as Figure 7 but for the eigenfunctions Sbi^ = kpahf / Bq (solid lines), Hbi^ = kpahf^ /Bq (long dashed lines), and 
= kpahj, /Bq (short dashed lines). 


the differential equation becomes singular at a point, and the situation we have in our calculations for low frequency range is 
essentially the same as the second order differential equation. 

This situation is largely different from that met by Lander et al. (2010), who found polar and axial Alfven modes for 
magnetized rotating stars with a purely toroidal background magnetic field. They suggested that pure Alfven modes of a non¬ 
rotating star or purely inertial modes of an unmagnetized rotating star may be replaced by hybrid magneto-inertial modes 
for magnetized and rotating stars (see also Mathis & Brye 2011, 2012), and that in the limit of MjT —>■ 0 or uAjCd —>■ 0 
the hybrid modes reduce to purely inertial modes, where M and T are the magnetic and kinetic energies of the equilibrium 
model, respectively. The method of calculation Lander et al (2010) use is a MHD simulation that follows the time development 
of linear oscillations around the equilibrium model and is different from the method we use in this paper. At this moment we 
do not understand why we find no magnetic modes for a purely toroidal background magnetic field. 


5 CONCLUSION 

In this paper, we have calculated non-axisymmetric oscillations of rotating stars magnetized with purely toroidal magnetic 
fields. Here we have used polytropic models of the indices n — 1, 1.5, and 3, and included the effects of the deformation caused 
by the toroidal magnetic fields. We have obtained discrete normal non-radial oscillation modes such as g, /, and p modes for 
non-rotating case. The frequency change due to the magnetic field for the modes scale with the square of typical stellar Alfven 
frequency uja = pB qKAtipq B?), and the proportional coefficients for the frequency changes are estimated by two different 
methods. For rotating stars, we have obtained rotational modes such as r and inertial modes. The frequency changes for the 
rotational modes also scale with the square of typical stellar Alfven frequency uja- From Tables 1 to 4, we find that high 
frequency modes such as / and p modes are susceptible to the stellar deformation, while the lower frequency modes such as 
g, r, and inertial modes are almost insensitive to the deformation. It may be important to note that in the present analysis 
we could not find jmax-independent magnetic modes, the existence of which are suggested by Lander et al (2010). The reason 
for the difference between the two calculations is not yet well understood. 

The present analysis is a part of our study of the oscillations of magnetized stars. Even for a purely poloidal magnetic 
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field, it is difficult to determine frequency spectra of non-axisymmetric modes and those of axisymmetric spheroidal modes. We 
note that the stability of a magnetic field configuration is another difficult problem. It is well known that a purely poloidal and 
purely toroidal magnetic fields are unstable and the energy of the field is dissipated quickly, that is, for several ten milliseconds 
(e.g., Goosens 1979; Kiuchi, Yoshida, & Shibata 2011; Laskey et al. 2011; Ciolfi & Rezzolla 2012), although stellar rotation 
may weaken the instability of a purely poloidal or purely toroidal magnetic fields (e.g.. Lander & Jones 2011a, 2011b). It is 
thus anticipated that a mixed poloidal and toroidal magnetic field configuration such as twisted-torus magnetic field (e.g., 
Braithwaite & Spruit 2004; Yoshida & Eriguchi 2006; Yoshida, Yoshida & Eriguchi 2006; Ciolfi et al. 2009) can be stable, that 
is, such a magnetic field configuration can last stably for a long time. If this is the case, it will be important to investigate 
the oscillation modes of stars threded by both toroidal and poloidal magnetic fields. As suggested by Colaiuda & Kokkotas 
(2012), however, the presence of a toroidal field component can significantly change the properties of the oscillation modes of 
magnetized neutron stars. In the presence of both poloidal and toroidal field components, toroidal and spheroidal modes are 
coupled, which inevitably breaks equatorial symmetry and antisymmetry of perturbations. 

As an important physical property inherent to cold neutron stars, we need to consider the effects of superfluidity and 
superconductivity of neutrons and protons on the oscillation modes (e.g., Andersson et al. 2009; Galmpedakis et al. 2011). It is 
believed that neutrons become a superfluid both in the inner crust and in the fluid core while protons can be superconducting 
in the core. For example, if the fluid core is a type I superconductor, magnetic fields will be expelled from the core region, 
because of the Meissner effect, and hence confined to the solid crust (e.g., Colaiuda et al. 2008; Sotani et al. 2008). In this 
case, we only have to consider oscillations of a magnetized crust so long as toroidal modes are concerned, and we have toroidal 
crust modes modified by a magnetic field, while spheroidal oscillations can be propagative both in the magnetic crust and in 
the non-magnetic fluid core. However, a recent analysis of the spectrum of timing noise for SGR 1806-20 and SGR 1900-1-14 
has suggested that the core region is a type II superconductor (Arras, Gumming & Thompson 2004). If this is the case, the 
fluid core can be threaded by a magnetic field and hence the frequency spectra of oscillation modes will be affected by the 
superconductivity in the core (e.g., Colaiuda et al. 2008; Sotani et al. 2008). 
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APPENDIX A: PULSATION EQUATIONS FOR THE SLOWLY ROTATING STAR WITH PURELY 
TOROIDAL MAGNETIG FIELDS 


To describe the master equations concisely, it is useful to introduce the following column vectors composed of the expansion 
coefficients for the perturbation quantities: the vectors S, H, , T, hS , and j/ji defined by 


{S)i=Si,, {T)i=T,,, {h% = hL {h^),=hh, 


Pi 

( 1 / 2 ). = 

pag 


(Al) 


where {X)j denotes the j-th component of the vector X and g = GM{a) jad is the gravitational acceleration. The perturbed 
continuity equation (1241) . and the a-, Q- and ^-components of the perturbed Euler equation (12311 . respectively, reduce to 


da 


Vg — 3 — a 


di9(a) 


da 


I — a 


d&ip) 

da 


Ao ^ S - Ucya + [Ao + 3i9(/3)Bo] H + 3mi?(/3)QoT, 


(A2) 


= [cicr^ {[1 -I- 277(0:)] I -f 2ri{l3)Ao} + aAll 8-^(2 + pcioj\ (I - Ao) . 

da'- ■'3\dlnay 


-|-(1 — aA — U)y 2 H —— ( 2 -I- 


d Inp 
din a 


pClUJA (I - Ao)y2 


—{277icicrfl[l -1-0-1- 77(0)]! -I- 2mciaQ[l3 + rj(/3}]Ao + Scia (3Bo}H 
—{2ciaCl[l + a + T7(a)]Co -I- 2ciGtl[j3 + rj{l3)]AoCo + I3 Qq}T 


, 1 - -2 , 

+ -^pciuJA \ rn 


da 


-h^ + 2{2 + 


d Inp 
d In a 


-Co 


dh^ r, /r, d Inp 

a— - 1-2 2 - 1 - ^ 


da 


din a 


- 2 

— {mu [1 + 0 - 1 - 77 ( 0 )] I + mu \P + vih)] Ao — 3/3 (2Ao + Bo)} 8 + aAp^ (2Ao + Bo) 8 

- —l^oy 2 — p-^Va ( 2 A 0 + Bo) 7/2 

cicT'' 

+ [(1 + 2 o) AqLo + 2/3 (AoAo + 3Bo) — 2muf3 (I + 6 A 0 )] H 
+ [-U (1 + 2o - 2/3) AoMi - 477/3 (AoAqMi + 3QoBi) + fiTn^gQo] T 


-2 


+ -mp— 2 + 


d In p 
din a 


A® + mp^h^ - = 0 , 

2 cr^ 


(A3) 


(A4) 
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-2 -2 

{v[l + a-\- J 7 (a)] AiK + i/ [/3 + ■q^P)] (AiAiK + SQiQoQi - 3Qi) - 3m,gQi} 5 - maAp^QiS + mp^VGQij/j 


+ [-u{l + 2 a- 2/3)AiMo - 4iy/3 (AiAiMo + 3QiBo) + 6 mPQi] H 
+ [(1 + 2a)AiLi + 2P (AiAi + 3Bi) - 2mi//3 (I + 6Ai)] T 




-nP^ 2 + 


d In p\ 
din ay 


-2 1-2 

AiKA® + p^^^AiMoh" + -mp‘LA (Ai - 21) = 0. 


The a-, Q- and ^-components of the perturbed induction equation (1261) . respectvely, lead 
hP = mS, 


(AS) 

(A 6 ) 


= maAI\Q — mVal^o ^?/2 + itt-H, 


(A7) 


lY = aAKS - VaKy^ - mT. (A 8 ) 

Here, 


_ d in M(a) 
din a ’ 


1 d Inp 
IT dlna’ 


■d(a) = 3a + 


da 
da ’ 


q[a) = a + a 


da 

da 


(A9) 


and (7 = aj(GMj is the frequency in the unit of the Kepler frequency, and v = 2 fl/(j. The quantities Qo, Qi, Co, Ci, 
K, Mo, Ml, Ao, Ai, R, Lo, Li, Ao, Ai, 6o, and Bi denote the matrices defined as follows: 

For even modes. 


{Qo)n = J!;+u (Qo)i+i,i = i*7+2, (Qi)ii=J7+i, (Qi)i,i+1 = J7+2, 


(Co)ii = -ih + 2)jr+ 1 , (Co)i+i.i = {ij + i)jr+ 2 , (Cl),- 


= IjJA+u (Ci ),,,+1 = -(Z,- + 3)jr+2, 


(K)ii 


Zi + 1 ’ 


(K)i,i+i 


jm 

Zj + 2 ’ 


(Mo),-,- — 


Zi + 1 


(Mo)j,j+i = 


Zi + 3 


Zi-+2 


i,-+2, 


(M 


iJit 


h + 2 

Ij + 1 




(Mi),-+i,,- 


Ij + 1 
Z,- + 2 


jm 

Ji,+ 2 , 


(Ao)ij — Ijih + l)i (Ai),-,- — (Z,- + 1)(Z,- + 2 ), 


^D^ di +2)(Z,- — 1) jm di + l)(^t +4) ,m 

—-;—r-;-A,+i, (.Kjj+i.j —- ; —"t;-A1+2, 


Zj + 1 


Z,- + 2 


Lo = I — mi/Ao Li = I — mul\i ^ = -(3QoQi — I), Ai = -(3QiQo — I), Bo — QoCi, Bi = QiCo 
where Z,- = \m\ +2j — 2 for j = 1, 2, 3, ....,ymax, and 


J m 
1 . = 


{Ij +m)(Z,- -m) 


1/2 


.(2Z,--1)(2Z,- + 1)_ 

For odd modes, 

(Qo),-,- = jr+i, (Qo)i,i+i = jr+ 2 , {Qi)n = JY+I, iQi)j+iJ = 


(AlO) 


(All) 


(Co),-,- = z,-jr+i, (Co),-.,-+i = -(Z,- + 3) Ji7+2, (Cl),-,- = -(Z,- + 2 ) jr+i, (Ci),-+i.,- = (z,- + i)jr+2, 


-'b+l 


J m 

i,-+2 


(K)ii =-(K),-+1 = 

^ z,- + l’ ^ Z,-+2 


(Mo),-,- — 


Ij + 2 


(Mo)i+i.i = ^^■^7+2, (Mi)ii = (Mi),-.,-+i = 


(Ao)ii — (Z,- + 1)(Z,- + 2), (Ai)j-j- — Ijilj + 1), 


(R)ii 


h ih A 3) 

Z,- + l 


(R)i,t+i 


h dj A 3) 

Z,-+2 


Lo = I — mvl\, 


‘0 


Li = I — mul\i 


A)= 2(3QoQi-I), 


Ai = -( 3 Q 1 Q 0 -I), 


Bo = QqCi, Bi = QiCo 


(A12) 
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where Ij = \m\ + 2j - 1 for j = 1, 2, 3,jmax- 

Eliminating the variables , and from equations (IA3ll - (IA5ll by using equations (IA 6 I) - (IA 8 I) . and eliminating H 

and T by using equations (lA4ll and (lA5l) . we may reduce equations (IA 2 | and (O to a set of coupled first-order linear ordinary 
differential equations for the functions y^ = S and j/j) which is formally written as: 


d 


da 


Vi 

?/2 



(A13) 


The surface boundary conditions are 


-V 1 +V 2 + 




= 0 , 


(A14) 


which means 


pV 


A p-b 


Stt 


|B| 


= 0 at the stellar surface, where AQ denotes the Lagrangian change of the quantity Q. 


The boundary conditions at the stellar center are the regularity conditions for the eigenfunctions j/j and 1/2 • 


APPENDIX B: SURFACE BOUNDARY CONDITION FOR THE FUNCTION ^>2 

Te determine the function 1^2 satisfying the differential equation IISI), we need the surface boundary conditions. Assuming 
the deviation of the surface r = Rs (R, 9) of the magnetized star from the surface r — R ot the non-magnetized star is small, 
we may write 

Rs{R,e) = R{l + 5<:{e)). (Bl) 

Since we have p{Ra, 6) = 0 and po{R) = 0 at the stellar surface, from equation ( 1101 ) we can obtain 

5C = [V'o(l) + V’ 2 (l)P 2 (cose)] . (B2) 

Using equation the gravitational potential ^{r,9) and its derivative d^{r,6)/dx inside the star are given by 

cE>(r, 9) = 4'o(r) -|- co - 2R'^uj\ [ci,o -b i)o{x) + V' 2 ( 2 :)P 2 (cos 6 ')] - [1 - P 2 (cos 6 ')] , (B3) 


d^{r,9) 

dx 


94'o(r) 

dx 


- 2 R^uj\ 


diJo , . dtp 2 . . . , 

(x) + -^{x)P 2 {cos 9 ) 


dx 





[1 - P2 (cos6 I)] , 


(B4) 


where we have set the constant C in equation (0 as C = co — 2P^a;ACi,o- Since p{Rb,9) — 0 and ’5o(Ps) ~ 4'o(P) + 
{d^ 0 /dr)r=RR5C, and (9To(r)/9r)r=fis = GM/R^ « {GM/R?){1 — 25C,) at the deformed surface, the gravitational potential 
^{Rs,9) and its derivative d^{Ra,9)/dx reduce to 

4?(Ps, = 4'o(P) + Co — 2P^a;ACi,o, (B5) 


a<E> 


a^o(p) 

dx 


- 2 R^uj\ 


1 ) + 


di )2 


('1 1 Po ("cos 01 


- [1 - P2(cos0)] - [V’o(l) -b V’2(l)P2(cOS0)] , 


(B 6 ) 


where we have used d^^/dx = GM/R? at the surface. On the other hand, the gravitational potential outside the star can be 
written as 


J, 1^0 od 2 2 1 ^ 1,0 , «:i ,2 „ , 

=- 2R uja —^ H-^P2(cos0) 

X I X X'^ J 


d^ _ KO 2 2 r«l,0 , o«:i,2 D , aO 
^ ^ + 2 ^? [^ + 3— P2(cos0)J , 

and at the stellar surface a: = Xs = 1 + (5C we have 

<I> = —Ko — 2R'^uj\ [ki,o + «:i,2P2(cos 0)] -b Kod(, 


(B7) 

(B 8 ) 

(B9) 


(yuj o o 

— = Ko + 2 R^cja [/^i,o + 3^1,2^2(0036')] — 2 ko^C^ (BIO) 

where kq, ki,o, and ki ^2 are arbitrary constants. If we require 4? and d^/dx inside and outside the star are continuous at the 
stellar surface, by comparing the zeroth-order terms, we find 


4'o(P) + Co = —Ko, Ko 


d'I>o(R) 


dx 


(Bll) 
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and hence we obtain 


^'o(^) — —Kq —- 


dx 


Co = 0. 


By comparing the perturbed terms we can obtain following relations: 
— Cl ,0 = —fvl ,0 + V'o(l), f^l,2 = 


(B12) 


(B13) 


dipo,.. ^dp 




(B14) 


From the relations, we note that the unknown constants ci,o and ki,o for the function tpo are determined uniquely by integrating 
equation m to the surface: 


dipo ^dp 


dLo f,-, , ^ dp 


On the other hand, for the function 'i/)2 we obtain by eliminating the constant «;i,2 


(B15) 




which gives the outer boundary condition for 'tp 2 at the surface. 


(B16) 


APPENDIX C: FREQUENCY CHANGES DUE TO THE MAGNETIC FIELDS 

Using the continuity equation (PI)l and the adiabatic relation (US, we may rewrite the Euler equation (1231) as 


- a^lil + 2 e) ^ + ae“Voe + a(^ • Voe)e,] = -Vox + e. —A 


Vo • $ + ^ • Vo ( 3 e + 


+ iaD 


_^^£j_Vo]n^_foVo_|) ^ ^ ^ ^ X B) X B' + (Vo X B') x B] 


(Cl) 


where x =p'!P- We write the eigenfunctions and eigenfrequency as follows (for a similar treatment, see, e.g., Saio 1981): 

^ = ^0 + ^2, (C2) 


X = Xo + X2, 


(C3) 


n = (To + (T 2 , 


(C4) 


where quantities with subscripts 0 and 2 denote quantities of order and ui\, respectively. The Coriolis term, B, is then 
written by 


where 


Bi°^ [^] = 2U sin [$] = 211 cos , B™ [^] = -211 (^sin 9^ + cos 


[^] = 2fl ( 2£ + ) sin Of, [^] = 211 ( 2e + 


de 


(2) rti _ 


Df[i] = -2^ 


da 


„ de\ . /„ sinScle 

2e + a_jsm0C +(2^+^^ 


sin 6 de 
cos 9 d9 

cos 9^^ 




Introducing equations (C2)-(C5) into equation (Cl) and grouping quantities of the same order in lja, we obtain 


Fip 


- o-Q^o = -VoXo + ea—AVo ■ ^o +'i<^oD^ 'Kol • 
for order oj^, and 


- (To^2 - 2(Toe^o “ f^oaCo Voe - (To(^o ' Voe)e,i - 2(Jo(72^o 


= -VoX2+e„^^A 
P 


^0 • ^2 + ^0 ■ ^0 ( 3 £ + 


+ iaoD^°^ [^2] + Kol + iaoB'-'^^ [^g] 


(0)r 


i(2)r 


_^ (^o-'^o Inp + Vo-^o) (y X B) X B + — [(V x B') x B + (V x B) x B'l 

47rp 47rp ( J 


(C5) 


(C6) 


(C7) 













Stars having purely toroidal magnetic field 17 


for order Note that Equation (C 6 ) for order describes the oscillation of the unmagnetized slowly rotating star. The 
fnnctions xo and X 2 , and B' are, in terms of and ^21 given by 

pFi 


Xo 


P 

pF: 


■ (Vo • ^0 + ^0 • Vo Inp - ^0 • SaA) 


A d €. \ 

X 2 = -- — -I Vo • ^2 + ^0 ■ ^0 ( 3e + ] + C 2 ' Vo Inp — ^2 ' > 


da J 




sin6»£i,„fc^oB™^ . 


sin 6 


(C 8 ) 


Mnltiplying Eqnation (C7) by the complex conjugate of the displacement vector and integrating over mass, we obtain the 
integral relation that includes no term related to ^ 2 ) given by 

pM pM pM pM 

- 2 ag / AMa - al {a^^V ot) ■ - al [(^o ' VoejeJ • ^SdM. - 2 aoa 2 / l^ol'dM. 

Jo Jo Jo Jo 

pM / c)f\ 

= xl ^0 • Vo [se + a—j AMa + iao [^ 0 ] • C^AMa + ia^ [^g] • ^JdM, 

1 pM 1 / ^ \ 1 pM -1 

+ ir - —^XO + ^0 • [(Vo X B) X B] • ^;dMa + — / - [(V x B') x B + (V x B) x B'l • (C9) 

47r Jo P \ pi 1 / Fit Jq p'- 

where AMa = p(a)a^ sinPdadPd((>. Substituting the 01 ( 4 -order eigenfunctions expanded like Equations (29)-(35) into Equation 
(C9) and taking 0-2 = we may obtain the integral expression for the coefficient E 2 , given by 

r 1 '^ 1 1 / 

E '2 = - / fi{a)ppa‘^Aa + ao / f 2 {a)pa‘^da +/ — f 3 {a)pa'^Aa + Ft / f 4 {a)pa'^Aa Wi, (CIO) 

.4ao Jo Jo 2cro Jg ci Jo J / 

where 
Wi = I 


+ i/iAo H + T+Ai t] pa'^Aa 
y [m + J/iS + |B|^ + |T|^) + S^Co T - T^AiK S + ffiAoMi T + TiAiMo b] 

d Inp 


pa^da, 


(Cll) 


dh 


fi (a) = S' < ma - mn + 2m ( 2 + 

' da ' 


+T'^ 


- 2 + 


d In p 
din a 

+Bi 

d Inp 


a"-C o 


m ( 2 + 


a^+2(2 + 


da 
d Inp 
d In a 


A" 


d In a 

+ 2mh^ — AoR A 


d In a 


AiK A® + 2AiMo a" + m(Ai - 2I)A^ 


2ayl 

'“3~ 


2 + 


d Inp 
din a 


Sfi\-Ao)S + 


2Vg 


2 + 


d Inp 
din a 


Sfi\-Ao)y2 


+2aAW {2Aq + Bo)S- 2VgH^ {2Ao + Bo) ^2 - 2maAT^QiS + 2mVGT^Qiy2 , 

/ 2 (a) = 5+ [p(d)l + y{P)Ao] S + a (^B+Ao B + T+Ai t) + ^ B^ {AqAo + 3Bo) B + PTfiAiAi + 3Bi) T 
-^/JB^Bo B + ^^B+ {2Ao +Bo)S + 3m^ (b+Qo T + T+Qi b) - (siQo T + TiQi s) , 


(C12) 


(C13) 


/3(a) = y2^ 


d'i9(a) diJ)/?) . 

a J M + a—fi^Ao 


da 


da 


S - Mfi)yfiBo H - 3my{P)y2^Qo T , 


(C14) 


/ 4 (a) = {[a + p(a)] I + [/3+ p(/ 3 )]A 3 } B - mBi {[a + p(a)] I + [/3 + p(/ 3 )]A 3 } B 

-2md (|B|2 + |T|2) - 4/3Ti (AiAiMg + 3QiBo) B 
-si {[d + p(d)]Co + 1/3 + p(/3)]A)Co} T - 2(d - P) BiAoMi T 
-2(d - /3) TiAi Mo B - 2m/3 Bi (I + 6 A 0 ) H - AP (A 3 A 0 M 1 + 3QoBi) T 

-2m/3 Ti (I + 6 A 1 ) T + [d + p(d)] TiAiK S + [/3 + p(/3)] Ti (AiAiK + SQiQgQi - 3Qi) S. (C15) 

Here, d = a/u)A, P = PISa, Hk = {GM/, and means the Hermitian conjugate of the complex matrix X. The 
magnetic perturbations , and A^ are given by 

= mS , 


(C16) 
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= maAI^Q ^ S — mVcI^Q ^ 1/2 + mH , 


(C17) 


= aAKS - VGKy 2 - mT. 


(C18) 


In the expression for E 2 given above, all the eigenfunctions, S, H, T, and 1/2 for unmagnetized stars even though the 
subscript “0” is not attached. 

When inertial modes are considered for the u;^-order eigensolution, for which lim ^ = kq with kq being a constant, 

—^0 a 2 

from Equation (CIO), we see that 


<^2 V2 -2 o , n 

— —^ as ^ 0, 

(Jo 


where 772 is a constant depending on the mode considered, which is given by 


V2 = - lim 


n-»o dfcH Wi 


pR pR 1 

/ }i{a)ppYda + 2 / —f3{<Ap' 

Jo Jo Cl 


•a'^da 


(C19) 


(C20) 


Equation f|C19[l implies that for the inertial mode, our expression for fJ 2 becomes inappropriate in the case of f2^ ^ oja- For the 
inertial mode, therefore, the condition ua f2^ ^ 1 is required for the expression for cr 2 to be applicable. Similar expressions 
to Eq. (|C19[) but for stars with general magnetic field distribution have been obtained by Morsink & Rezania (2002). 
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